ROM-based Quantum Computation: Experimental Explorations using 
Nuclear Magnètic Resonance, and Future Prospects 
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ROM-based quantum computation (QC) is an alternative to oracle-based QC. R has the advan- 
tages of being less "magical", and being more suited to implementing space-efíicient computation 
(i.e. computation using the minimum number of writable qubits). Here we consider a number 
of small (one and two-qubit) quantum algorithms illustrating different aspects of ROM-based QC. 
They are: (a) a one-qubit algorithm to solve the Deutsch problem; (b) a one-qubit binary multi- 
plication algorithm; (c) a two-qubit controlled binary multiplication algorithm; and (d) a two-qubit 
ROM-based version of the Deutsch- Jozsa algorithm. For each algorithm we present experimental 
verification using NMR ensemble QC. The average fidelities for the implementation were in the 
ranges 0.9 - 0.97 for the one-qubit algorithms, and 0.84 - 0.94 for the two-qubit algorithms. We 
conclude with a discussion of future prospects for ROM-based quantum computation. We propose 
a four-qubit algorithm, using Grover's iterate, for solving a miniature "real-world" problem relating 
to the lengths of paths in a network. 

PACS numbers: 03.65.Yz, 03.75.Fi, 42.50.Lc, 03.65.Ta 
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I. INTRODUCTION 

The current excitement in, and perhaps even the 
existence of, the field of quantum computation |ÏJ] is 
due to the demonstration that quantum computers can 
solve problems in fewer steps than classical computers 
l@j H> Q ü' Hi' improvement is rigorously established 
for the Deutsch- Josza algorithm [[| and for Grover's 
search algorithm Q , while Shor's factorization algorithm 
Q uses exponentially fewer steps than any known classi- 
cal algorithm. 

It is interesting that, of the above quantum algorithms, 
those that are provably faster (a) are not exponentially 
faster, and (b) make use of an oracle. An oracle is a 
"black-box" that defines a function 
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(1.1) 



Here Z^v is the natural numbers modulo N, that is, 
{0, 1, . . . , N — 1}. The oracle Of acts on a n-qubit string 
\k), and an m-qubit string \l) as follows: 



O f \k)\l) = \k)\l®f(k)), 



(1.2) 



where ® represents bit-wise addition modulo 2. Noto 
that we are defining an oracle so that it can be applied 
to classical bit strings as well as to qubit strings. 

Although the concept of an oracle is very useful in the 
context of complexity theory, they are, as their name 
suggests, somewhat "magical" in their operation. Thus 
they may hide a great deal of computational complex- 
ity in one step, and for this reason can be considered 
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"unrealistic" j7j. In a quantum context it has been sug- 
gested that counting oracle calls may be a poor way to 
study the power of algorithms ||. Finally, it seems to 
us that oracle-based computing is best for studying time 
efficiency, rather than space eflicieney 

All of these factors suggest that it is worth explor- 
ing an alternative basis for computation. In this paper 
we explore quantum computation based on ROM (Read- 
Only Memory). In an earlier paper two of us and co- 
workers showed that a ROM-based quantum computer 
is more space-efíicient than a ROM-based classical com- 
puter. Here space eflicieney is defined in terms of the 
number of writable qubits required. In particular, one 
writable qubit is sufficient to compute any binary func- 
tion of an arbitrary number of ROM bits, whereas two 
writable bits are needed to achieve the same. Also, for 
a particular one-bit function (multiplication of all the 
ROM bits) evidence was found to support the conjecture 
that one qubit can solve the problem in polynomial time, 
whereas three bits are required for the same. 

These results indicate that ROM-based computation is 
ideal for demonstrating space- (and possibly time-) effi- 
ciency on small scale quantum computers. Of course, at 
the moment small scale quantum computers are all we 
have experimentally. For example, in ion traps the num- 
ber of qubits that can be coherently controlled is at most 
four pfj| , and in Nuclear Magnètic Resonance (NMR) 
experiments on ensembles of molècules, the number is 
at most seven [jn]. In this paper we explore the space- 
efíicient quantum algorithms in Ref. [^) , as well as other 
ROM-based quantum algorithms, in an NMR context. 

The structure of this paper is as follows. In Sec. II we 
review ROM-based computation as defined in Ref. || . In 
Sec. III we present the simplest space-efíicient one-qubit 
algorithm (which solves Deutsch's problem). Sec. IV 
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covers the one-qubit ROM-multiplication algorithm of 
Rcf. ||. In Sec. V we present a two-qubit version of this, 
the controlled-ROM-multiplication, which is also prov- 
ably more space efíicient than any classical algorithm 
(which would rcquire three bits). In Sec. VI we exploro 
the Deutsch-Josza algorithm using ROM rather than an 
oracle. In each of the sections III-VI we present ex- 
perimental results following the theory, and we discuss 
these results in Sec. VII. We conclude in Sec. VIII with a 
discussion of future prospects for ROM-based quantum 
computation, and in particular we propose a four-qubit 
demonstration of quantum computing solving a "realis- 
tic" problem (i.e. a problem that can be related to the 
real world and that would require more than a second of 
human thought to solve). 



II. ROM-BASED COMPUTATION 

We consider quantum computation using qubits, the 
numbcr of which rcmains fixed throughout the computa- 
tion. The computer evolves by the operation of gates, 
which implement a unitary operation on one or more 
qubits simultaneously. It has been shown |Ï2] ] that a 
single two-qubit gate, such as the controlled-NOT gate, 
supplemented by all one-qubit gates, is sufficient to per- 
forin all possible quantum computations in this model. 
Unitary gates are of course reversible. This means that 
in principle the computation can be carried out without 
dissipation of information and hcnce without energy cost 

§0- 

To make a fair comparison with unitary quantum com- 
putation, we must consider reversible classical compu- 
tation. As is well known, universal reversible classical 
computation is not possible with just one-bit and two-bit 
gates. Rather, a three-bit gate such as the Toffoli gate 
or Fredkin gate is required p] , The measurement of 
the state of the qubits (in the computational basis) takes 
place only at the end of the computation. Similarly, ini- 
tialization (setting a bit to a fiducial state such as |0)) is 
allowed only at the beginning of the computation. These 
stipulations are necessary to keep the computation non- 
dissipative. 

Before proceeding, let us establish some notation. We 
will write the n-bit (or qubit) representation of a num- 
ber a; £ Z 2 ™ as \x). This is equivalent to the notation 
|x) = |x„_i)|x„_ 2 ) • • ■ |xi>|ar ), where x = J2 p x p 2P - In a 
'circuit' diagram, the most significant bit (MSB), |a; n _i), 
will appear at the bottom of the diagram, and the lcast 
significant bit (LSB), \xo) at the top. 

We used this notation already in Eq. (|l.2| ) to spec- 
ify the action of an or acle which implements the func- 
tion / defined in Eq. (1.1). In ROM-based computa- 
tion, the function / that is the subject of the compu- 
tation is implemented not by an oracle, but by its vàl- 
ues {/(fc) : k G Zjv} being stored in rea d-only memory. 
Specifically, for / as defined in Eq. (1.1), N x m ROM 
bits are required to store the function. For the simple 



case m = 1 (a binary function), we require N bits which 
could be allocated as /o,/i, ■ • • , /jv-i ; where = f(k). 
These ROM bits are not counted in the size of the com- 
puter. That is to say, the size of the computer is taken 
to be the number of additional (non-ROM) (qu)bits. 

To capture the essence of read-only memory, we impose 
the following constraints: 

1. The ROM bits {/& : k} can be prepared only in a 
classical state. 

2. For any gate involving the writable qubits, any sin- 
gle ROM bit, fk for some fc, may act as an addi- 
tional control bit. 

3. No other gates involve the ROM bits. 

These three conditions together imply that the ROM bits 
will always remain in the same state. In finite state au- 
tomata models, space-bounded computation can be dis- 
cussed using Turing machincs with two tapes, one of 
which is read-only J?J. This is clearly very similar to 
the present idea of individually-accessed ROM bits. The 
necessity for placing a constraint on the number of ROM 
bits that can act as simultaneous control bits was dis- 
cussed in Ref. || . 

The restriction to single-bit ROM access leads to a 
simplification in the representation of ROM in circuit di- 
agrams of reversible computation. Rather than explicitly 
using "wires" to represent the ROM-bits we will simply 
leave a space at the top of the diagram, and write in 
which ROM bit (if any) is acting as the extra control bit 
for that gate. This suggests an alternative way to concep- 
tualize the replacement of the oracle by ROM. An oracle 
is like an all-knowing person who refuses to divulge in- 
formation except whcn askcd a question in a certain way. 
ROM is like a committee of people who each have one 
bit of information but who refuse to communicate with 
one another except by acting individually upon a device. 
In this way problems in ROM-based quantum computa- 
tion can be seen to have some similarities to problems in 
quantum communication such as in Refs. ÍÏJm fL7|. 



III. ONE QUBIT SOLUTION TO THE 
DEUTSCH PROBLEM 



A. Theory 

The smallest quantum computer is obviously one 
qubit. It turns out that this, plus additional ROM bits 
rather than an oracle, is sufficient to solve the Deutsch 
problem Q for any n. The Deutsch problem can be 
phrased in the following way. Given a function of the 
form ([□]) with N = 2™ > 4 and m = 1, find a true 
statement from the following list: 

(A) / is not balanced. 

(B) / is not constant. 
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A constant function / is one for which ^ fc f(k) = or 
N; that is, for which f(k) = V k or f(k) = 1 V k. 
A balanced function / is one for which ^2 k f(k) = N/2. 
Clearly one of (A) and (B) must be true, and they both 
may be true in which casc cither can be chosen. 

Deutsch and Josza found a quantum algorithm that 
solved this problem using n+1 qubits and two oracle calls 
j|. By replacing the oracle with 2™ ROM bits, we are 
able to solve the problem with a single qubit and with one 
control from each ROM bit. If we were concerned with 
time-efficiency, the exponential number of "ROM calls" 
may seem a problem. However here we are concerned 
only with spacc-efhcicncy. 

The one-qubit algorithm to solve this problem is very 
simple: 



fo h 



Ín-i 



(3.1) 



The computer is prepared in the fiducial state |0). Each 
ROM bit, /fe G {f , fi,--- , /iv-i}i m turn controls (indi- 
cated by the vertical line) a rotation on the qubit with 
unitary operator [^] . That is, the gate is implemented 
if and only if /fe = 1 . Here we are using the notat ion 



[6] a = exp[-i(6/2)* a }, 



(3.2) 



where a a are the usual 2x2 Pauli matrices, with a G 
{x, y, z}. For the Standard representation of these matri- 
ces, the basis states are 



|0) = 



II) = 



(3.3) 



In Eq. ( fr.l| ) , the measurement is represented symbolically 
by an eye: £- , and yields the result x, a single bit. That 
this is a classical piece of information is represented by 
the double, rather than single, wire. 

If the function is constant, then either it never leaves 
the state |0), or it is rotated by N x (2tï/N) = 2tt around 
the y axis, returning it to the state |0). If the function 
is balanced, it is rotated by (N/2) x (2tï/N) — ir around 
the y axis, putting it into the state |1). If it is ncithcr 
balanced nor constant it will end up in a superposition 
of |0) and |1), so a measurement will yield either result. 
This computation clearly solves the Deutsch problem. If 
the measured state x of the computer is 0, the answer 
returned is (B). If the measured state is 1, the answer 
returned is (A). 

To show the superiority of a space-bounded quantum 
computer over a space-bounded classical computer we 
simply have to prové that a one-bit classical computer 
cannot solve the D J problem. Consider the simplest case, 
where n — 2, so that / maps {0, 1, 2, 3} to {0, 1}. Since 
the only possible one-bit gate is a NOT gate [N], which 
obeys [N] 2 = 1, the only one bit operation for this prob- 
lem is 



where each pk G {0, 1}. Acting on the initial state 0, this 
computes the functional "^Pkfk modulo 2. It is triv- 
ial to prové that this functional does not distinguish be- 
tween balanced and constant functions for any choice of 

P0,Pl,P2,P3- 



B. Method 

The sample used for all of the following experimental 
demonstrations was a 0.1M solution of heavy chloroform, 
"C^Clg, dissolved in d6-acetone, CD 3 COCD 3 (for lock- 
ing purposes). Chlorine isotopes 35 Cl and 37 C1 have large 
quadrupole moments (I = 3/2), resulting in extremely 
short relaxation times when covalently bonded, on the 
order of 10/i,s. This has the effect of mas king scalar cou- 
pling between chlorine and other nuclei | ^9[ . Thus, the 
chloroform molecule is effectively a two-spin system, pro- 
ton and carbon-13, with 1=1/2 for both spins. 

All spectra were obtained using a Bruker DRX-500 
spectrometer, for which the magnitude of H$ was approx- 
imately 11. 6T. The resonance freqüències of the proton 
peaks were Un = 500.137849MHz and u c = 125.77754749 
MHz. The scalar coupling was measured to be J = 
(214.8 ± 0.5)Hz. Clearly, J < - i/ c |, so that the two 
spins can be resonantly excited independently. There is 
more than 1kHz separation to the solvent lines, which 
thus played no part in the experiment. All experiments 
were performed at a temperature of (298 ± 0.1)K. The 
measured vàlues for the longitudinal T± and transverse T 2 * 
(including field inhomogeneity effects) relaxation times 
were ïi(H) = (9.7 ± 0.2)s, 7\(C) = (11.0 ± 0.2)s, 
T 2 *(H) = (6.4±0.3)s, and T 2 *(C) = (0.2 ± 0.01)s. The 
maximum pulse program time was approximately 20ms, 
significantly less than all of the above vàlues. 

For the one qubit algorithms, the H nucleus was used 
as it had a far narrower linewidth. The initial state is 
the thermal equilibrium state, which has a small excess 
spin in the longitudinal direction (spin up) . This pseudo- 
pure state |2(| has observable signal proportional to that 
of state |0), as desired. A one-qubit gate can be imple- 
mented by an appropriately phased transverse magnètic 
field pulse (or short sequence of pulses), rotating at the 
resonant (radio) frequency of the nucleus. If a particular 
gate is ROM-controlled then it is implemented only when 
the value of the controlling ROM bit is one. Following 
the complete pulse sequence, a tt/2 transverse pulse is 
used to shift longitudinal spin into the transverse plane, 
where its precession will induce a signal in the RF coils 
(the read-out). A positive spectrum indicates an excess 
of spin-up populations before the read-out pulse was ap- 
plied. The presence of the 13 C nucleus, with almost equal 
population spin up and spin down, causes a frequency 
splitting of J/2. Thus the observed spectra for the two 
logical states are of the form 



[N] 



/0P0+/1P1+/2P2+/3P3 



(3.4) 



10) 



-rr (3.5) 
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The fidelity of the transformation is calculating by di- 
viding the area under the spectrum by that which would 
have arisen from a perfect transformation of the thermal 
signal. Since the final readout is equivalent to the aver- 
age of the results of projective measurements in the a z 
basis of each member of the ensemble, the area ratio R 
can be considered to be due to a mixture of the correct 
result (with probability F) and the incorrect result (with 
probability 1 — F), namely R = F - (1 - F). Thus the 
fidelity is calculated asF=(iï+l)/2. 



C. Results 

The Deutsch problem has a deterministic output if one 
adds the promise that the function / is either balanced 
or constant. In this case output (A) indicates that / is 
constant and (B) that it is balanced. With TV = 4, this 
means that in effect there are only three different pulse 
sequences arising from the algorithm in Eq. (3.1): that 
in which all vàlues of / are zero, that in which two are 
one, and that in which all four are one. 

The results of these three different pulse sequences are 
shown in Fig. 0. We see that the results agree well with 
the theory. The first case consists of doing nothing, so 
its fidelity is one, by definition. The fidelity of the other 
cases is calculating as described above. The average fi- 
delity (taking into account that there are six possible 
ways in which the function can be balanced) is F = 0.9. 



on the fact that if 2" mod 3 = 2 then 2™" 1 mod 3 = 1 
and vice versa.) 

In this section we consider another one-qubit algo- 
rithm, which is also impossible on a one-bit computer, 
and which is conjectured to be more time efficient 
than any two-bit algorithm. It also solves a more nat- 
ural problem than the DJ problem, namely to find the 
product of N ROM bits u±, ■ ■ ■ u^. The quantum algo- 
rithm derived in Ref. [^) requires exactly ./V 2 ROM-calls 
for N a power of two, and 0(N 2 ) otherwise. The required 
number of ROM-calls r for a two-bit classical computer 
was found by numerical search to be r — 1,3,5,9 for 
N = 1,2,3,4. It is conjectured that r(N) is given by 
the recursion relation r(N) = r(N — 1) + 2L W / 2 J, and 
there is an obvious classical algorithm requiring exactly 
this many ROM-calls. This formula is clearly asymptot- 
ically exponential in n, but is actually smaller than the 
TV 2 ROM-calls in the quantum algorithm for N = 2, 4, 
and 8. 

In the experiment we only implemented the quantum 
algorithm for N — 2 and N = 4. The one qubit algorithm 
which determines u±xu2 can be constructed as follows. 



Ul U2 U\ U 2 



(4.1) 



In an abuse of our notation, we will indicate the above 
algorithm as 



0000 
F=1 



0011 
F=0.89 



1111 
F=0.86 



FIG. 1: Spectra for H nucleus showing the implementation of 
our one-qubit solution to the Deutsch problem. The vàlues 
of the four ROM bits are shown above each spectrum. The 
fidelities (F) are also shown. 



Ui X íí 2 



(4.2) 



Similarly, a gate effecting the transformation , 
conditional on u\ x u%, is 



Ui U2 Ui U2 Ui X IÍ2 



(4.3) 



These operations can be combined to construct an algo- 
rithm which determines the answer a — ui x u 2 x u 3 x u±, 
viz. 



IV. ONE QUBIT MULTIPLICATION 



A. Theory 



Ui XU2 U3X U4 Ui XU2 U3X 1*4 

1°)- [- f ] — W — [i] — M — £= 



(4.4) 



The above results demonstrate that a one-qubit quan- 
tum computer can solve problems that no one-bit classi- 
cal computer can. By adding one more classical bit, the 
problem can be solved (this is always true, as shown in 
Ref. P]). Moreover, there is a two-bit algorithm that is 
just as time efficient as our quantum algorithm above. (It 
uses the two bits to tally the fkS modulo 3, and is based 



B. Results 

The results shown in Figs. and ^ were obtained by the 
method outlined above. Again we see good agreement 
with theory. The average fidelity was F — 0.97 in the 
case of multiplying two ROM bits, and F — 0.92 for 
multiplying four ROM bits. 
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00 



F=1 



01 

F=0.98 



10 

F=0.99 



11 

F=0.91 



FIG. 2: Spectra for the H nucleus for the one-qubit algorithm 
for multiplying two ROM bits (shown as M1112). A small sys- 
tematic error is evident in the dispersive features seen in the 
last case. Other details are as in Fig. 1. 



the value of the second writable bit. The modified al- 
gorithm is identical to the one in the preceding section, 
except that all of the gates (which act on the first bit 
\xo)) are controlled by the second bit x\ as well as by all 
of the ROM bits. That this cannot be done on a two-bit 
classical computer follows from the proof by Toffoli |ÏS|| 
that multi-controlled-NOTs can not be built from single 
controlled-NOTs without the use of an auxiliary writable 
bit. 

For the case n = 2, the circuit is 

Ml U 2 Ui U 2 

l°Hf];rMir[-fL-WirS= ^xu.xu^. (5.1) 

\x\) • • • • — £-= x\ 



F = l 
0000 




F=0 . 99 
0001 




F=0 .89 
0010 




Here the sòlid circles on the wire for the second qubit 
indicate that it acts as a control qubit for the relevant 
gate. 


F=0 . 92 

0011 

1 


J 


F=0 . 97 
0100 




F=0 . 92 

0101 

1 - 




B. Method 


F=0 . 91 

0110 


J 


F=0 . 90 
0111 




F=0 . 91 

1000 




As well as selective RF pulses tuned to the H and C 
nuclei, the two-qubit algorithm requires an interaction 
between the nuclei. This occurs simply by leaving time 


F=0 . 91 
1001 




F=0 . 90 
1010 




F=0 . 98 

1011 




between the (negligibly short) pulses for the spin-spin 
coupling Hamiltonian 

H = hJafa^/i (5.2) 



F=0 . 90 






F=0 . 92 




F=0 . 91 




1100 






1101 




1110 





to act. These periods of free evolution are usually of 
duration 1/4 J or 1/2 J, and are simply denoted by this 
time. For example, 



F=0 ._8_5 
1111 



1 

2J 



exp[-iH(l/2J)/h] 



1 



(5.3) 



FIG. 3: Spectra for the H nucleus for the one-qubit algorithm 
for multiplying four ROM bits (shown as U1M2W3U4). Small 
systematic errors are evident in the dispersive features seen 
in most cases. Other details are as in Fig. 1. 



V. TWO QUBIT 
CONTROLLED-MULTIPLICATION 

A. Theory 

It is simple to generalize the above single-qubit algo- 
rithm (with a space advantage of one bit compared to 
classical computation) to a two-qubit algorithm also with 
a space advantage of one bit. This is done by requir- 
ing the calculation of u\ x u 2 x ■ ■ ■ u n x x%, where x\ is 



The two-qubit algorithm also requires a pure initial 
state. A pseudo-pure state 1 00} of both spins (H and C) 
up can be prepared from the thermal equilibrium state 
by a sequence of RF pulses, free evolution, and gradient 
pulses. The last of these effectively removes the trans- 
verse spin of the sample. That is, it diagonalizes the state 
matrix into the logical basis. We use the pulse sequence 
of Cory et al. [ pï[ , but change the [ir/6] y pulses for both 
spins into a [— 7r/6L pulse. This is to ensure that the 
signal is that for the pseudo-pure state |00), rather than 
the negative signal, which corresponds to a state matrix 
oc I — |00)(00|, where / is the 4x4 identity matrix. In 
theory, this pseudo-pure state preparation procedure re- 
sults in a signal reduction by 3/8 compared to the original 
thermal equilibrium state. 

The readout is done identically to the one qubit case. 
The correspondence between the logical states and the 
observed spectra is more complicated. The single reso- 
nance peak for each spin is potentially split into a dou- 
blet, at ui ± 7t J . With the NMR convention of frequeney 
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increasing from right to left, the spectral shapes for the 
four logical states are as follows. 



State Hydrogen Carbon 

|00) — L -1 
|01) 
|10) 



T 



(5.4) 



111) 



As in the one qubit cases, for calculating the fidehty 
we are interested only in the occupation of the logical 
states, since our computation is meant to be determin- 
istic. We integrate under the spectra at the four fre- 
qüències, and divide by 0.375 of the area of the original 
thermal state. The 0.375 is due to the theoretical sig- 
nal loss in the pseudo-pure state preparation described 
above. This procedure gives a number linearly related to 
the occupation probabilities poo, í>io> Poi, an d Pn- For 
example, the area ratio i?p under the left Hydrogen peak 
should satisfy 



Ri =PiQ-Pu· 



(5.5) 



In addition, the probabilities should sum to unity. Thus 
we have five equations in four unknowns, which we solve 
by a least-squares method to yield the probabilities. If 
the desired outcome is |11), for example, then the fidclity 
equals p\\. 

The above algorithm requires four controlled gates. 
These can be constructed from resonant pulses on the 
two nuclei, plus periods of free evolution of duration 1/2 J 
or 1/4 J. Specifically, 



i 4 ''> 



ï\x [ 2 ] ^ T ~l [siy [ i\ x 



(5.6) 



(5.7) 



and 



L 3 Jn 



(5.8) 



they should be since the control (Carbon) qubit being 
set to zero means that Xq = x\ x u\ x u 2 = 0. The last 
four spectra are repeats of the first four, but with the 
control (Carbon) qubit initially rotated from |0) to |1). 
This state, |10), was prepared with fidehty 0.89, as shown 
for the case where the ROM bits are 00. As expected, all 
spectra but the last also have xo = 0, and the last shows 

Xq = 1. 



Hydrogen 



Carbon 




where n is an axis defined by n = (x + y)/ 



|00>,00 
F=0.93 

|00>,01 
F=0.83 

|00>,10 
F=0.92 

|00>,11 
F=0.83 

|10>,00 
F=0.89 

|10>,01 
F=0.81 

|10>,10 
F=0.85 

|10>,11 
F=0.76 



FIG. 4: Spectra for the H and C nuclei for the two-qubit 
algorithm for multiplying two ROM bits, controlled by the 
first (C) bit. The initial pseudopure states of the two bits 
are shown with the spectra, along with the ROM bits (shown 
as u\U2). Small systematic errors are evident in most cases. 
Other details are as in Fig. 1. 



C. Results 



The pseudo-pure state |00) was produced with fidelity 
of 0.93. This appears as the first line of Fig. [|, which is a 
running of the controllcd-multiplication algorithm when 
the ROM bits are 00 (i.e. nothing is done). The results 
for the other possible ROM vàlues appear in the next 
three spectra. All four of these spectra are the same, as 



The raw fidelities, calculated as discussed above, are 
shown on the figure. Let us scale out the fidelity of the 
initial state preparation (that is, divide the raw fidelities 
by 0.93 for the first four spectra and 0.89 for the second 
four), and take the average. Then we get a mean fidelity 
for the implementation of the algorithm of F — 0.94. 
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VI. TWO QUBIT ROM-BASED 
DEUTSCH-JOSZA ALGORITHM 



A. Theory 



case N = 4 (n — 2), the algorithm is 



We saw in Sec. [II that the Dcutsch [|| problem 
can be solved on a ROM-based computer with a single 
qubit. This algorithm was quite unlike that proposed by 
Deutsch and Deutsch and Josza ||. In this section 
we investigate the implementation of the Deutsch- Josza 
algorithm on a ROM-based computer. This requires at 
least two bits to solve the Deutsch problem. Our motiva- 
tions here thus do not include space efficiency. Instead, 
they are as follows. 

First, as noted in the introduction, ROM-based com- 
putation seems more realistic, so it is interesting to see 
how it can be applied to an apparently quintessentially 
oracular algorithm. 

Second, there is a question of interpretation of past 
experiments. Again as noted in t he i ntroduction, an or- 
acle should be definable [see Eq. ( |l.2|) 1 by its action on a 
classical computer. This is necessary in order not to give 
an unfair space advantage (of m qubits) to a quantum 
computer. This requirement is met in the original the- 
oretical proposals of Deutsch and Josza || and Grover 
[^J. However, it is not met in proposals such as that in 
the "refined" Deutsch- Josza algorithm of Ref. |^2| , implc- 
mented in Ref. J23| . That is because in this algorithm the 
oracle directly produces phase shifts, whi ch h ave no clas- 
sical analogue. The requirement of Eq. (1.2) would also 
rule out the oracles implemented in other NMR experi- 
ments [|§ H, ü ( but not t0 th °se in Refs. §7], §§, f§). 
Our analysis here will show that these experiments can 
be very easily reinterpreted in terms of ROM calls rather 
than oracle-calls. 

Third, there is a question of how quantum the Deutsch- 
Josza algorithm is. The use of a non-classical oracle al- 
lows the Deutsch- Josza algorithm to be implemented us- 
ing one fewer qubit (n rather than n + 1). The same 
number (n) of qubits are required for the ROM-based 
implementation. For the minimal case n = 2, it was 
shown in Ref. |2^] that with the non-classical oracle, the 
Deutsch- Josza algorithm does not utilizc entanglement. 
On this basis, the authors claim that it therefore "solves 
the Deutsch problem in a classical way." Leaving aside 
qüestions as to the meaning of "classical" in this context, 
we show that in a ROM-based implementation, entangle- 
ment necessarily occurs. This suggests that the so-called 
classicality noted in Ref. [^2| is due to the unrealistic 
nature of the oracle they use. 

In the ROM-based implementation of the Deutsch- 
Josza algorithm, the ROM bits are the same as in Sec. III, 
namely the binary vàlues /o, • • • , }n-i of the function /, 
which is either balanced or constant. For the minimal 



|o>- 
|o>- 



-[f] F 



/oo 

1 




/oi 

1 




/io 

1 




/ll 
1 


000 




001 




010 




011 

















X . 



Xi 

Here the four distinct two-qubit gates, <f> a p change the 
sign of the logical state |a/3), leaving the other three un- 
altered. Mathematically, the operation of these gates can 
be expressed as 



4>a(}\l8) = (1 - 28 ai 8ps)\iS) 



(6.2) 



The result x\Xq = 00 indicates a constant function; any 
other result indicates a balanced function. 

The four ROM-controlled gates together have exactly 
the same effect as the non-classical oracle introduced in 
Ref. |^2). This is an example of how any non-classical 
oracle has a ROM analogue. The application of an odd 
number of these gates creates an entangled state, so the 
intermediate states of the quantum computer are entan- 
gled. It is only because the number of times the gates are 
applied is even (because the function is promised to be 
balanced or constant) that the state at the end of the four 
ROM-controlled gates is not entangled. Entanglement is 
required in all cases of this algorithm except when the 
function is identically zero. 



B. Method 

We use the following realization of the two-qubit phase 
ïates: 









0Q/3 




1 




2J 



(6.3) 



C. Results 

The res ults are shown in Fig. ^. Comparing with the 
table (5T), we see that the result 00 is obtained only 
for the case of function vàlues 0000 and 1111, as ex- 
pected. Results for the vàlues 1100, 1010 and 0110 were 
not obtained. If the algorithm implementation were per- 
fect then the state after the four controlled phase gates 
in these cases would be the same as for 0011, 0101, and 
1001 shown in Fig. 0. There is also no reason from the 
pulse sequence to expect the fidelities to be very different 
for those cases implemented. We can thus take the cases 
shown as being representative, and use thcm to calcu- 
late an average fidelity over the eight possible functions 
/. Scaling away the fidelity of 0.93 for the pseudo-pure 
state preparation (as in Sec. V), we obtain F = 0.84. 
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Carbon 



Hydrogen 



0000 
F=0 .86 



0011 
F=0 . 75 



0101 
F=0 .81 



1001 
F=0 .76 



1111 
F=0 . 75 



FIG. 5: Spectra for the H and C nuclei for the two-qubit 
ROM-based Deutsch-Josza algorithm. The ROM bits are 
shown as /oo/oi/io/n- Note that the order of the Carbon 
and Hyd rogen spectra are reversed as compared to Fig. ^] and 
Eq. ( |5.4|) . Small systematic errors are evident in all cases. 
Other details are as in Fig. 1. 



VII. DISCUSSION OF EXPERIMENTAL 
RESULTS 



The experimental results shown clearly verify the 
ROM-based quantum algorithms discussed in this paper. 
The average hdelities of the algorithms were not umisu- 
ally good for NMR experiments. They were 0.9, 0.97 and 
0.92 for the one qubit algorithms, and 0.94 and 0.84 for 
the two-qubit algorithms. These hdelities are dehnitely 
correlated with the length of the pulse sequence required. 
However, it is interesting that the lowest fidelity (0.84) 
was obtained for the the ROM-based Deutsch-Josza algo- 
rithm, which was not much longer than the other 2-qubit 
algorithm, but which differed from it in that used entan- 
gling operations. 

The largest source of error was probably spatial and 
temporal variation in the intensity and phase of the RF 
pulses. The spatial variation of field strength is a di- 
rect consequence of limitations imposed by the structure 
of the coils, being small Hclmholtz coils. Cummins and 
Jones Q provide a good discussion of the errors incurred 
in NMR computing and explain how the systematic er- 
rors due to Hi held inhomogeneities can be greatly re- 
duced. We did not attempt to apply these techniques. 



VIII. FUTURE PROSPECTS: A "REALISTIC" 
PROBLEM 

In this paper we have presented a number of ROM- 
based quantum algorithms for one- and two- qubit pro- 
cessors. One of these (one qubit multiplication) was 
derived in Ref. |J; the rest are new. They include 
a one-qubit algorithm solving the Deutsch problem, a 
two-qubit controllcd-multiplication algorithm, and a two- 
qubit ROM-based Deutsch-Josza algorithm solving the 
Deutsch problem. For all algorithms we have also pre- 
sented experimental verification, using NMR enscmble 
quantum computing. 

All bar one of the above algorithms demonstrated 
space efficiency, in that a classical computer would re- 
quire an extra processor bit to solve the problems. The 
exception is the ROM-based Deutsch-Josza algorithm. 
We believe that future prospects for ROM-based quan- 
tum computation lie more in the direction of this last 
example. There are two reasons. First, it follows from 
the results of Ref. that the maximum space efficiency 
offered by ROM-based quantum computation is one bit, 
which could not be significant in computations of a use- 
ful scale. Second, in showing how an oracular algorithm 
can be implemcnted on a ROM-based computer, the ex- 
ample of Sec. VI illustrates how time-efficient quantum 
computing could be implemented realistically. 

In the remainder of this section we will explore a future 
prospect for ROM-based quantum computation along 
these lines. We take as our basis not the Deutsch-Josza 
algorithm, but the other famous oracle-based algorithm, 
due to Grover We will show how the oracle in this 
algorithm can also be implemented using ROM-calls. We 
find a specific implementation with two properties of in- 
terest. First, it is experimentally feasible in the short 
term, requiring only four qubits. Second, it solves a prob- 
lem that can be related to a read-world situation (albeit 
a miniaturized one), and that would require more than a 
second of cerebral processing time to solve. 

The problem we consider relates to the lengths of paths 
between two vertices in a network (a set of vertices con- 
nected by edges of diffcring lengths). Some problems 
of this nature, such as hnding the longest such path, are 
known to be NP-complete JjJ . This is a class of problems 
that are almost certainly exponentially hard to solve, and 
are thus of great practical interest. 

Consider the network bclow 



U 



B 



\ 



E 



(8.1) 



There are four vertices, labeled B, L, U, and E (for Be- 
gin, Lowcr, Uppcr, and End), linked by edges. These 
could represent cities and roads respectively. We are in- 
terested in the length of paths from B to E, as indicated 
by the direction of the arrows in Eq. (8.1). Assuming no 



back-tracking, there are four possible paths: BLE, BUE, 
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BLUE, and BULE, to which we assign the numbers from 
to 3. Each path p has a length L{p) associated with 
it, equal to the sum of the lengths l(e) of each edge e 
of the path. Six edges must be distinguished, as the 
edges UL and LU could well have different associated 
lengths. This is because "length" could represent some 
generalized cost, such as the time a traveler would have 
to wait for a lift. We discretize the problem by assuming 
that for all e, l(e) is either zero or one. Thus for all p, 
L(p)€ {0,1,2,3}. 

The obvious question a traveler would like to ask is, 
what is the shortest path? Unfortunately, this is not the 
sort of question that Grover's search algorithm will an- 
swer straight away. Rather, what it can answer is qües- 
tions like, what is the path of length 1? If there is ex- 
actly one path of length 1, Grover's algorithm will find 
it. If there are none or two paths of length 1, Grover's 
algorithm will return one of the four paths at random. 
If there are three paths of length 1, Grover's algorithm 
will actually return the only path that is not of length 1! 
Thus, what Grover's algorithm really does in our case, is 
to return a number p which means that 

path p, or none of the paths, is the odd-one- 
out with respect to having the length L. 

Here the odd-one-out is the only one having, or the only 
one lacking, a property. 

With a little thought it is apparent that actually we are 
not limited to making a demand about a specific length 
L. Rather, we can demand information about a set § of 
lengths L. Since the total length L e {0,1,2,3}, there 
are seven such sets, 

Si = {0}, (8.2) 
§2 = {1}, (8.3) 



§3 


= {o,i}, 


(8.4) 


§4 


= {2}, 


(8.5) 




= {0,2}, 


(8.6) 


Se 


= {1,2}, 


(8.7) 


§ 7 


= {0,1,2}. 


(8.8) 



There are other nontrivial subsets of {0, 1, 2, 3}, but they 
are the complements of the above seven sets, so they 
would lead to demands already covered by the above 
seven sets. Specifically, the seven demands we could 
make on our computer are, with j € {1, • • • , 7}, 

What is a path p such that p, or none of the 
paths, is the odd-one-out with respect to hav- 
ing a length L(p) £ §j? 

In a ROM-based computation, there are six ROM bits to 
encode the lengths 

Ibl, Ibu, Ilxj, Ile, and ?ue- (8-9) 

and three to encode one of the seven sets §j , 

jo, ji, and j 2 , (8.10) 

where these are the bits in the binary representation of 
j. The vàlues of the nine ROM bits thus code for 448 
different instances of the general problem. 

Note that the binary representation J2J1J0 of j is re- 
lated to the set §.,■ as follows: 

^ = {L : JL = 1} (8.11) 

Using this we can construct a solution to the above prob- 
lem by the following circuit: 



^UE jo jl h ^BL 



«UE 



|0>— 

|o>— 

I0MH]- 



C/bl 







000 




001 




010 






t/UB 
















^BL 









The first pair of Hadamard gates jïj puts the 1 0000) 
state into the state 



(1/2)5>>|00>, 

p=0 



(8.13) 



where the lower two qubits, encoding the path, are in 
a superposition of all four possible paths. The next 
six gates, conditioned on the six ROM bits storing the 



17, 



UE 



HO) 
HO) 

H77H 
H77H 



.12) 



'ou 



HtíH^= po 



r 



lengths of the edges, use the upper pair of bits to count 
the length of each path. For example, 



U B j J \ P )\L) = \p)\L + X [BLr P }). 



.14) 



Here x[H] ^ s the characteristic function, equal to 1 if its 
argument H is true, and zero otherwise, while "BL □ p" 
means "BL is an edge in path number p" . These func- 
tions are 'hard-wired' in the gate construction, as they 
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depend only of the geometry of the netw ork i n Eq. (3.1), 
not the lengths. The addition in Eq. flS.14 ) is defined 
modulo 4, which is reversible on 2 bits, so that this is a 
wcll-dchncd gatc. After all six such gates have acted, the 
state of the processor is 



(1/2) j>)|L(p)>. 



(8.15) 



p=0 



Upon this superposition of all paths, and their associ- 
ated lengths, we now act the sign change that is at the 
heart of Grover's algorithm. The three gates controllcd 
by j = 4j2 + 2ji + jo produce the overall phase shift 



\p)\L)^(í-2 X [Le§ J })\p)\L). 



(8.16) 



This c hange s the sign of the components of the superpo- 
sition ( 8.15 ) for which the length is in the specified set Sj. 
After the application of the inverse of the six controllcd 
gates Ubl • • • C^ue, the processor is in the state 



(l/2)£{l-2x[i(p)GS i ]}|p)|0). 



.17) 



p=0 



Past experimental implementations of Grover's algo- 
rithm [^D, ^5[ [2(| have relied upon a non-classical or- 
acle that takes the proces sor d irectly from a state like 
Eq. ( $.13 ) to one like Eq. flB.17 ). In that case, the upper 
pair of qubits in Eq. fl8.12| ) are of course superfluous. For 
these experiments, where there is no real problem being 
solved, such an oracle seems reasonable enough. How- 
ever, in the context of the present problem, which relates 
•'real-world" situation ■ 



.1 



an or- 



to a "real-world" situation — the network 
acle like this would indeed be magical. The point of the 
above analysis is to show that, with the help of just two 
additional qubits, the required oracle can nevertheless be 
implemented in a realistic manner, using ROM-calls. 

A Grover iterate is complete with the application of 
Hadamard gates to the lower pair of qubits, and a phase 
change to the 1 00) state p3|. In this case, because the 
number of paths is four, a single Grover iterate suffices. 
The final step of the algorithm is to apply the Hadamard 
gates again. After this, the state of the lower pair of 
qubits of the processor is \pj), where p — pj is the uniquc 
solution of L(p) 6 Sj or L(p) G Sj, where is the com- 
plement of Sj . If neither of these equations have a uniquc 
solution, then the final state is a superposition of all pos- 
sible paths, as in Eq. ( S.13 ). Thus it is apparent that 
this algorithm does indeed fulfill a demand of the form 
above. 

The above algorithm is certainly not space-efficient. 
From the results of Ref. || it follows that even a clas- 
sical computer could solve this problem with a two-bit 
processor (although probably in more steps). Nor do we 
claim that it is time-efficient. A classical four-bit proces- 
sor may well be able to solve the problem in fewer steps. 



However it would be interesting to determine what the 
effect would be if one stipulated that the three demand- 
specifying ROM-bits (the bits of j) only be used onec, 
as in the above quantum algorithm. A similar 'once- 
only' constraint on information access was considered in 
Ref. @ 

A lack of both time and space efficiency for this partic- 
ular algorithm would not render it worthless. Considcr 
the one-qubit ROM-multiplication algorithm in Sec. IV. 
The specific instances of that algorithm we discussed 
and experimentally implemented were not time-efficient. 
However, when sealed up to larger numbers of ROM 
bits, the algorithm was (we conjectured) time-efficient 
compared to the minimal (two-bit) classical processor 
needed to solve this problem. Similarly, for large prob- 
lems, a suitable generalization of this ROM-based Grover 
algorithm would, we hope, become quadratically time- 
efficient compared to any classical algorithm. 

If Grover's algorithm cannot be applied to "real-world" 
problems in this way, then it is of very limited utility. 
Arguments pointing to the generality of its quadratic 
speed up for large problems have been made by two 
groups. First, Brassard, H0yer and Tapp showed 
how Grover's algorithm can work even if the number of 
"marked" elements is unknown. Second, Cerf, Grover 
and Williams claimed that Grover's algorithm can 
make use of the structure of a large scale problem in the 
same way as classical search algorithms, and thus main- 
tain a quadratic speed up. However, serious doubts on 
this matter have also been expressed p4| , so the question 
remains open. 

Experimentally implementing the above four-qubit 
algorithm is well beyond the scope of this study. We 
have not even compiled it into the appropriate "machine 
language" (NMR pulse sequences). However, this could 
be done along the same lines as the other algorithms 
presented he re, b y first decomposing the four-qubit 
gates in Eq. ( $.12 ) into sequences of one and two-qubit 
gates. We hope that the challenge of experimental im- 
plementation on this, or some similar algorithm, is taken 
up. Although it would still be only a toy calculation, it 
would be a significant step on the way towards full-scale 
calculations of difficult problems pertaining to real-world 
situations. 
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